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Abstract-  This  paper  describes  a  magnetic 
anomaly  guidance  system  that,  with  support  from  the 
Office  of  Naval  Research,  is  being  developed  for  fully 
autonomous  detection,  localization  and  classification  of 
ferrous  mines  in  Very  Shallow  Water/Surf  Zone 
(VSW/SZ)  environments. 

The  magnetic  guidance  system’s  hardware 
configurations  and  magnetic  target  signature  processing 
methods  specifically  have  been  developed  for  autonomous 
guidance  of  small,  maneuverable  sensing  platforms,  such 
as  Underwater  Bottom  Vehicles,  (UBVs),  to  magnetic 
targets.  The  prototype  magnetic  anomaly  guidance  is  a 
three-axis  magnetometer-gradiometer  comprised  by  an 
array  of  triaxial  fluxgate  magnetometers  and  triaxial 
accelerometers.  This  design  allows  development  of 
magnetic  gradient  tensor  contraction  scalars  that  provide 
the  basis  for  a  unique  method  for  magnetic  guidance  of 
UBVs  to  underwater  and  buried  mines.  The  triaxial 
accelerometer  array  can  be  used  for  mitigation  of  the 
target-detection-range-reducing  effects  of  large  changes 
in  sensor  platform  orientation. 

Field  test  results  indicate  that  a  simple  magnetic 
guidance  algorithm  based  on  comparison  of  the 
magnitudes  of  contractions  of  gradient  tensor  subsets 
measured  by  “primary  guidance  axes”  of  the  sensor  array 
will  provide  a  robust  and  practical  method  for  fully 
autonomous  guidance  of  high-  mobility  UBVs  to  magnetic 
targets.  Target  detection  range  can  be  enhanced  by  using 
rotationally  invariant  tensor  contractions  for  robust 
proximity  sensing  and  “pseudo-robust”  scalar  total  field 
quantities  for  “homing  in  on”  magnetic  targets.  Sensor 
system  stability  and  DLC  performance  also  would  be 
improved  by  the  use  of  better,  temperature-compensated 
magnetometers. 

The  results  indicate  that  the  magnetic  anomaly 
guidance  system  and  its  method  for  target  localization  are 
particularly  appropriate  for  homing  in  on  magnetic 
targets  from  high-mobility  sensing  platforms  such  as 
small  UBVs  or  Navy  Divers. 


I.  INTRODUCTION 

With  support  from  the  Office  of  Naval  Research  (ONR), 
NAVSEA -Coastal  Systems  Station  (CSS)  has  participated  in 
the  development  of  magnetic  sensing  technologies  that 
potentially  will  enhance  the  effectiveness  of  autonomous 
underwater  vehicles  (AUVs)  for  mine  countermeasures 
(MCM)  in  littoral  Very  Shallow  Water/Surf  Zone  (VSW/SZ) 


environments.  In  particular,  CSS  has  developed  an  innovative 
magnetic  anomaly  sensing  approach  for  detection, 
localization  and  classification  (DLC)  of  magnetic  targets  such 
as  underwater  and  buried  mines  [1,2,3],  The  approach 
processes  target  magnetic  anomaly  signatures  and  develops 
magnetic  scalar  quantities  that  form  the  basis  of  a  unique 
magnetic  “Scalar  Triangulation  and  Ranging  (STAR)”[2] 
method  for  DLC  of  mines  using  high  mobility  sensing 
platforms  such  as  AUVs  or  Navy  divers.  In  particular  the 
magnetic  STAR  method  can  be  adapted  for  guidance  of 
autonomous  underwater  bottom  vehicles  (UBVs)  to 
underwater  mines  in  challenging  VSW/SZ  environments. 

This  paper  presents  field  test  data  from  a  prototype 
magnetic  anomaly  guidance  system  for  UBVs.  The  data 
indicate  that  magnetic  sensor  designs  based  on  the  Scalar 
Triangulation  and  Ranging  concept  also  would  provide  very 
effective  mine  DLC  systems  for  other  high-mobility  sensing 
“platforms”  such  as  Navy  Explosive  Ordnance  Disposal 
(EOD)  Divers. 

Generally  speaking,  the  physical  basis  for  the  use  of 
magnetic  anomaly  sensors  for  DLC  of  underwater  and  buried 
mines  is  the  magnetic  fields  BA  that  emanate  from 
magnetically  polarized  mines.  The  BA  fields  constitute  dc 
magnetic  anomalies  within  the  Earth's  background  field.  The 
magnetic  anomaly  fields  or  “target  signatures”  contain 
implicit  information  that  relates  to  the  location  and  magnetic 
strength  of  the  mines.  The  target  signatures  essentially  are  not 
affected  by  water  turbulence,  turbidity,  bubble  clouds,  multi- 
path  propagation  effects  or  even  intervening  layers  of  mud  or 
sand  (for  buried  mines)  that  can  limit  the  performance  of 
other  sensing  technologies  in  VSW/SZ  environments. 

However,  there  are  several  challenges  to  practical  use  of 
magnetic  sensors  for  DLC  of  mines  from  small,  high- 
mobility  sensing  platforms  such  as  UBVs  or  Navy  EOD 
Divers: 

•  UBV  operation  in  turbulent  VSW/SZ  environments 
typically  will  result  in  large  changes  in  sensor  platform 
orientation  that  can  result  in  reduction  of  DLC  range  due 
to  non-target-related  changes  in  vector  magnetic  field 
measurements.  Mobile  sensor  applications  generally 
either  require  devices  that  measure  magnetic  field 
gradients,  i.e.,  “gradiometers”  or  devices  such  as  scalar 
magnetometers  that  are  relatively  insensitive  to  platform 
orientation.  However,  conventional  gradiometry  usually 
requires  fairly  uniform  platform  motion  and  scalar 
magnetometry  does  provide  an  efficient  target  DLC 
modality  for  UBVs. 
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•  The  small  physical  size,  power  and  computational 
budgets  of  UBVs  largely  precludes  the  use  of  more 
conventional  tensor  gradiometry  and  scalar 
magnetometry  approaches  for  target  localization  and 
magnetic  anomaly  guidance. 

•  The  UBV's  self  magnetic  signature  due  to  electric  drive 
motors  and  other  electrical  subsystems,  will,  unless 
mitigated,  tend  to  overwhelm  the  very  small  target 
signatures,  BA  and  seriously  degrade  the  practical  DLC 
capabilities  of  any  UBV-mounted  magnetic  sensor 
system.  Consequently,  the  UBV  signature  should  be 
reduced  as  much  as  possible  and  residual  signature 
effects  should  reduced  by  judicious  application  of 
adaptive  signature  compensation  techniques. 

The  Scalar  Triangulation  and  Ranging-  STAR  concept  for 
DLC  of  magnetic  targets  was  specifically  designed  to  meet 
the  challenges  presented  by  small,  high-mobility  UBV-type 
sensing  platforms.  Initial  results  from  a  prototype  sensor 
demonstrated  the  validity  of  the  STAR  method  for  three- 
dimensional  DLC  of  magnetic  targets  [3].  This  paper  presents 
additional  test  results  that  more  explicitly  demonstrate  that 
the  gradient-contraction-based  STAR  concept  is  particularly 
appropriate  for  adaptation  into  a  relatively  simple  sensor 
system  for  fully  autonomous  guidance  of  UBVs  to  magnetic 
targets. 

On  continuation,  this  paper  will  present  the  following: 

1.  A  brief  technical  background  of  magnetic  anomaly 
sensing  theory  and  techniques. 

2.  Technical  details  of  the  magnetic  anomaly  guidance 
system  and  its  Scalar  Triangulation  method.. 

3.  Results  and  discussion  of  field  tests  of  the  sensor  system 
and  its  magnetic  guidance  parameters. 

II.  TECHNICAL  BACKGROUND 

In  general,  detection,  localization  and  classification  (DLC) 
of  magnetic  object  relates  to: 

•  The  existence  of  vector  magnetic  fields  BA  that  emanate 
from  the  objects’  net  magnetization  M.  The  BA  fields 
constitute  “magnetic  anomalies”  within  the  Earth's 
background  field. 

•  The  relatively  large  and  locally  nearly  constant  magnetic 
field  of  Earth,  BE.  Since  the  relatively  small  target 
anomaly  field  is  convolved  with  BE  ,  the  Earth's  field 
also  complicates  the  process  of  DLC  of  magnetic  objects 
using  mobile  sensing  platforms. 

At  distances,  r,  greater  than  about  three  times  a  magnetic 
object’s  longest  dimension,  the  “far-field"  BA  of  a 
magnetically  polarized  object  is  given  by  the  well-known 
magnetostatic  dipole  equation. 

BA(r,M)  =  (|i/4it)[  3(M»r)r/r5  -  M/r3]  (2.1) 

In  SI  units,  BA  (r)  [Tesla,  T]  is  the  vector  magnetic  induction 
at  a  point  defined  by  a  position  vector  r  [meters,  m]  relative 
to  a  vector  magnetic  dipole  source  M  with  units  of  Amperes 
x  meter  squared.  Am2.  The  parameter  u  [Tm/A]  is  the 


magnetic  permeability  which  for  non-magnetic  media,  is 
approximately  4ji  x  10’7  Tm/A,  and  it  =  3.14159... 

Direct  measurements  of  BA  are  complicated  by  the 
presence  of  the  relatively  very  large  Earth  field  BE  .  BE  has  a 
nearly  constant  magnitude  BE  of  about  («)  50,000  nT  (1  nT  = 
10"9  T)  at  the  Earth’s  surface.  However,  (2.1)  indicates  that 
dipole  fields  decrease  with  the  inverse  cube  of  distance  to  the 
objects,  i.e.,  BA  oc  r.  In  particular,  the  BA  fields  of  mine-like 
targets  drop  off  to  nT  and  sub-nT  levels  at  distances  of  a  few 
meters  from  the  targets.  Therefore,  BA  «  BE  except  for  field 
points  that  are  measured  very  close  to  the  target.  Also,  direct 
measurements  of  magnetic  field  actually  involve 
measurement  of  a  total  field  BT  that  is  the  vector  sum  of  BA 
and  Be,  i.e.,  BT  =  BE  +  BA.  The  magnitude  of  BT  in  Cartesian 
(XYZ)  coordinates  is 

Bt  =  (Btx“  +  Bty  2  +  Btz2)0'5  (2.2) 

Bt  is  a  rotationally  invariant  scalar  that  is  not  affected  by 
changes  in  mobile  sensor  platform  orientation.  Unfortunately 
however,  BT  doesn't  provide  a  robust  basis  for  target 
localization  since  in  some  regions  of  target  space  BT  may 
decrease  as  the  target  is  approached.  This  non-robust 
behavior  arises  because  locally  the  BA  field  lines  curve  360° 
through  space  while  BE  is  essentially  constant  in  direction. 
Consequently,  in  some  regions,  components  of  the  BA  fields 
point  in  the  same  direction  as  BE  and  increase  the  magnitude 
of  Bt,  in  other  regions  the  BE  and  BA  fields  point  in  opposite 
directions  thereby  decreasing  the  magnitude  of  BE. 

While  Bt  by  itself  is  non-robust,  the  Scalar  Triangulation 
and  Ranging  method  does  develop  a  BT  -based  “pseudo- 
robust”  [2]  quantity  that  can  be  used  for  longer  range 
magnetic  anomaly  detection  and  guidance  than  the  robust 
gradient  (G)-based  guidance  parameters  that  are  described  in 
the  following  paragraphs.. 

The  magnetic  gradient  tensor  and  magnetic  gradiometry. 

The  gradient  of  BE  is  quite  small,  typically  V  BE  «  0.02 
nT/meter.  Therefore,  the  effects  of  BE  can  be  greatly  reduced 
by  using  sensors  that  measure  the  gradient  of  BT  (i.e., 
“gradiometers”)  since  V  BT  »  V  BA.  In  magnetic 
gradiometry  the  measurements  of  gradient  components  of  BA 
provide  the  basis  for  target  localization  and  classification. 

The  gradient,  G  =VBA,  ( |  VBA(r)  |  k  r4  )  of  the  vector 
field  of  (2.1)  is  a  second  rank  tensor  whose  matrix  elements 
are  given  by  the  following  expression  [4]. 

(VB)jj  =Gii  =  3Bi/Sr7-  = 

=  -3((.i/4jrj[M.r(5r,r;  -  r%j)  -  r2(r,M j  +  r,Mi)]  r1  (2.3a) 

Eq.  (2.3)  constitutes  a  set  of  equations  that  relate  measured 
gradient  components  with  the  unknown  components  of  target 
position  r  and  dipole  moment  M.  In  (2.3a),  the  r,  ,  M,  terms 
respectively  refer  to  the  X,Y,Z  Cartesian  components  of  r 
and  M.  (r,  =  X,  r2  =  Y,  r3  =  Z  etc.).  Magnetic  gradiometry 
basically  involves  measurements  of  gradient  components 
combined  with  rather  complex  methods  for  gradient  tensor 
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“inversion”  to  determine  components  of  r  and  M  [4].  The 
complete  gradient,  Gt  matrix  is: 

Gxx  Gxy  Gxz  0Bx/0x  5Bx/5y  0Bx/0z 

G  s  Gyx  Gyy  Gyz  =  dBJdx  dBJdy  cBJcz 
Gzx  GzyGzz  0Bz/0x  dBJdy  dBJdz 

ABX/Ax  ABX/Ay  ABX/Az 
«  ABy/Ax  ABy/Ay  ABy/Az 

ABZ/Ax  ABZ/Ay  ABZ/Az  (2.3b) 

The  finite  difference  elements,  ABX/Ax,  ABX/Ay,  ABX/Az,  etc, 
represent  the  fact  that,  in  practice,  the  gradient  components, 
dBJdx,  dBJdy,  0Bx/5z,  etc.)  are  determined  by  measuring 
vector  field  components  at  spatially  points  along  the  XYZ 
axes,  subtracting  the  field  values  at  one  point  from  the  values 
at  another  and  dividing  the  resultant  differential  field  values 
by  the  separations  AX,  AY.AZ  between  the  measuring  points. 

For  conventional  gradiometry  the  AX,  AY.AZ  are  much 
smaller  than  the  distance  to  target,  r.  However,  the 
operational  paradigm  of  mine  DLC  using  UBVs  requires  the 
UBV  to  approach  to  within  such  close  proximity  to  the  mines 
that  the  approximations  G  oc  r  1  and/or  AX,  AY.AZ  «  r  may 
not  be  valid. 

Some  simplification  of  the  gradient  tensor  measurement 
process  can  result  from  Maxwell's  Equations  for  the 
divergence  and  curl  of  B,  namely,  V»B  =  0  and  V  x  B  =  0. 
Then  the  sum  of  diagonal  matrix  elements  is  zero  (i.e., 
0Bx/0x+SBy/5y  +  dBJdz  =0)  and  “across-diagonal”  matrix 
elements  are  equal  (that  is  dBJdy  =  dBJdx,  dBJdz  =  dBJdx, 
and  dBJdz  =  dBJdy)  and  measurement  of  five  independent 
tensor  components  will  completely  determine  the  gradient 
tensor.  The  present  work  uses  the  aforementioned  symmetry 
properties  to  construct  the  complete  gradient  tensor  GT  and 
the  corresponding  total  gradient  contraction,  CT  . 

However,  even  with  the  reduced  number  of  independent 
matrix  components,  the  localization  of  magnetic  targets 
through  inversion  of  the  gradient  equations  (2.3)  typically 
requires  operational  constraints  on  sensor  platform  motion 
(e.g.,  a  requirement  for  constant  velocity)  in  addition  to 
mathematical  ambiguities  and  requirements  that  G  oc  f 4  and 
AX,  AY.AZ  »  r,  that  may  make  the  gradient  tensor  inversion 
approach  more  complex  than  is  necessary  for  simple 
guidance  of  UBV-type  sensing  platforms  or  Navy  Divers  to 
magnetic  targets. 

A  simplified  approach  to  magnetic  anomaly  guidance  can  be 
based  on  the  Scalar  Triangulation  and  Ranging  (STAR) 
method.  The  STAR  method  develops  target  localization 
information  from  sums-of-squares  or  “contractions”  of  the 
gradient  tensor  components  [1-3].  In  Cartesian  coordinates,  the 
mathematical  operation  of  tensor  contraction  of  (2.3)  is  given 
by  the  following  expression: 

Ct2=  2  (G,j)2  =  Gxx2  +  Gxy2  +  . . . 

10 

*  (ABx/AX)2  +  (ABx/AY)2 +...  (2.4) 


The  contraction,  Ct  of  GT,  the  full,  nine  component  gradient 
tensor,  ,  is  a  rotationally  invariant  and  robust  quantity.  Here, 
rotational  invariance  implies  that  a  quantity  doesn't  change  as 
the  sensor  platform  is  rotated  relative  to  the  target. 
Consequently,  by  itself,  a  fully  rotational  invariant  quantity 
such  as  CT  does  not  provide  explicit  directional  control 
parameters  that  will  immediately  indicate  when  an  UBV  turns 
away  from  the  target  direction.  However,  as  will  be  shown  in 
the  EXPERIMENTAL  RESULTS  section,  partial,  3- 
component  subsets  of  the  magnetic  anomaly  gradient  tensor 
can  indeed  provide  a  robust  basis  for  guidance  of  UBVs  or 
Navy  Divers  to  magnetic  targets. 

III.  TECHNICAL  APPROACH 

The  need  for  a  compact  and  practical  magnetic  anomaly 
guidance  (or  homing  )  system  for  small,  high-mobility  UBVs 
has  led  to  a  simplified  adaptation  of  the  STAR  approach  [2,3] 
for  magnetic  anomaly  localization.  The  simplified  guidance 
approach  is  based  on  the  following: 

(A)  A  magnetometer-gradiometer  sensor  system,  shown  in 
Figs.  1  and  2,  that  combines  signals  from  a  three  axis, 
array  of  commercially  available  triaxial 
magnetometer/accelerometer  (TMA)  field  sensing 
elements  [3]  .  Each  TMA  contains  a  triaxial  fluxgate,  a 
triaxial  accelerometer,  A/D  converter,  microprocessor, 
and  serial  data  RS232  interface.  The  magnetometers 
measure  vector  components  of  B-field  at  each  point  of 
the  array.  The  accelerometers  measure  sensor  platform 
orientation.  Right  and  Left  “Primary  Guidance  Axes”,  R 
and  L  in  Fig.  2,  are  symmetrically  disposed  about  a 
central  symmetry  axis  S.  For  UBV  guidance,  magnetic 
signals  from  the  primary  guidance  axes  provide  left-right 
directional  control  for  “homing-in-on”  magnetic  targets. 

(B)  A  processor  uses  the  TMA  data  to  develop  magnetic 
gradient  tensor  and  scalar  total  field  quantities  that 
implicitly  contain  target  DLC  information.  In  particular, 
partial  three-component  gradients  (Gj)  are  calculated  for 
each  axis  (i.e.,  each  pair  of  TMAs).  The  partial  gradients 
are  used  to  construct  the  full  gradient  tensor,  GT.  Primary 
magnetic  anomaly  guidance  parameters  Cr,Cl  and  Ct  are 
developed  through  tensor  contraction  of  the  G,  and  Gt 
gradient  tensor  matrix  elements.  “Pseudo-robust”  sums 
of  total  fields  (2.2)  BTR,  BTl  and  BTT  are  also  developed 
for  the  sensor  axes  to  provide  longer  range  detection  and 
guidance. 

Specific  magnetic  guidance  parameters  are  developed  in 
the  following  manner.  Each  TMA  measures  the  X.Y.Z 
components  of  the  local  magnetic  field  BT.  The  BTX,  BTY,  and 
BTz  data  are  processed  to  develop  gradient  and  total  field 
parameters.  Gradient  components  are  calculated  for  each 
sensor  axis  in  the  vehicle  frame,  namely  ABXl/d,.  ABYl/d,,  and 
ABzi/di  (where  i  denotes  the  respective  R.L.B  and  S  axes  and 
dj  is  the  axis  length).  The  gradient  contraction  C,  is  calculated 
for  each  of  the  three  axes  using  (2.4).  For  example, 

CR  =  [(ABXR/dR)2  +  (ABYR/dR)2  +  (ABZR/AdR)2f5  defines  a 
partial  gradient  contraction  of  the  three  gradient  components 
that  are  measured  by  TMAs  2  &  3. along  the  sensor  system’s 
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R-axis.  The  “total”  gradient  contraction  Ct  is  given  by  the 
sum  of  squares  of  the  full,  nine  component  gradient  tensor, 
Gt-  Six  components  of  Gt  directly  are  obtained  from  direct 
field  measurements  and  three  components  are  obtained  from 
the  symmetry  properties  of  G.  Similarly,  total  field  values 


BTi,  BT2  and  BT3  are  calculated  for  each  of  the  three  TMAs 
using  (2.2).  Then  total  field  sums  are  calculated  for  each  axis. 
The  resulting  pseudo-robust  guidance  parameters  are  BTR  = 
Bt2  Bt3  ,  Bjl  =BTi  +Bt2  and  Btt  =  Bjj  +Bt2  +  Bt3. 


Fig.  1.  Prototype  gradiometer  and  magnetic  guidance  system.  Triaxial  Magnetometer-Accelerometers 
(TMAs)  measure  magnetic  field  and  acceleration  at  each  point  of  the  array.  The  Right  and  Left  (R  and 
L)  axes  are  the  “Primary  Guidance  Axes”  of  the  gradiometer. 
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Fig.  2.  Schematic  of  the  geometrical  configuration  of  the  prototype  magnetic  anomaly  guidance  system. 
Three  TMAs,  1,2,&3,  are  located  at  the  vertices  of  an  isosceles  triangle.  The  plane  of  the  triangle  is 
perpendicular  to  the  horizontal,  XY  plane  of  the  circular  sensor  support  platform.  The  Right  (R)  and  Left 
(L)  Primary  Guidance  Axes  respectively  are  comprised  TMAs  2&3  and  TMAs  1&3.  The  R  and  L  sides 
are  each  33.8  cm  long.  A  27.0  cm  symmetry  axis  (S)  goes  from  TMA-3  to  the  midpoint  (p)  of  Base  Axis 
(B).  TMAs  1  &2  comprise  the  40.3  cm  B-axis.  The  processor  develops  gradient  contraction  quantities  that 
are  used  in  the  Scalar  Triangulation  and  Ranging  (STAR)  method  of  magnetic  target  localization.  For  two 
dimensional  localization  of  mines  using  this  simple  three  axis  array,  the  Primary  Guidance  Axis  (R  or  L) 
that  is  closest  to  target  develops  the  largest  gradient  contraction  signal,  Cr  or  Cl.  Homing  in  on  a  target 
simply  involves  maneuvering  the  UBV  under  control  of  the  Primary  Guidance  Axis  signals  to  maintain  Cr 
=  Cl  and  verifying  increasing  values  of  the  Total  gradient  contraction,  Ct,  as  the  UBV  approaches  the 
target. 
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Extensive  field  tests  were  performed  to  verify  that  the 
magnetic  STAR  concept  will  provide  a  robust  basis  for  fully 
autonomous  guidance  of  UBVs  to  magnetic  targets.  The  test 
range  and  procedures  are  indicated  in  Fig.  3.  A  magnetic 
dipole  target  was  configured  to  approximate  the  anomaly 
field  be  developed  by  a  0.027  m3  ferrous  object  in  the  Earth's 
magnetic  field  [5,6],  Data  runs  were  performed  with  the 
target  near  the  plane  of  the  gradiometer  support  structure  or  1 
m  below  the  base  axis  of  the  sensor.  Consequently,  the  data 
represent  the  target  signatures  that  an  UBV  would  see  for  a 
mine  located  on  or  below  the  seafloor.  Previous  results  [3] 
indicated  that  the  sensor  system  was  susceptible  to  “gradient 
drift''.  In  order  to  quantify  a  putative  temperature  dependence 
for  the  sensor  drift  problem  TMA  temperature  data  were 
taken  and  processed  simultaneously  along  with  the  magnetic 
anomaly  data.. 

IV.  RESULTS  AND  DISCUSSION 

The  data  curves  in  Figs.  4-9  represent  gradient  contraction 
and  total  field  sums  from  the  Left  and  Right  “Primary 
Guidance  Axes”  (i.e..  Cl,  Cr,  Btl,  Btr,),  total  gradient 
contraction,  (Ct)  and  the  total  field  sum  from  all  three  TMAs 
(Btt).  Figures  4-6  show  the  variations  in  gradient  contraction 
and  total  field  values  that  were  measured  by  the  gradiometer/ 
guidance  system  as  it  was  moved  on  tangential  paths  with 
headings  of  135°,  180°  and  270°  and  with  closest  point  of 
approach  (CPA)  of  1.5  m  as  shown  in  Parts  (b)  and  (c)  of  Fig. 
3.  For  these,  and  other  headings  in  Fig.  3(b),  at  the  CPA  the 
gradiometer’ s  Left  (L)  axis  was  about  20  cm  closer  to  the 
target  than  the  Right  (R)  axis. 


All  magnetic-data  curves  are  shown  with  the  initial  “offset” 
value  at  the  beginning  of  the  data  run  subtracted  from 
subsequent  data  points. 

The  data  in  Figs.  4-6  show  that  within  the  gradiometer 
detection  range  of  about  2  m,  CT  >  CL  >CR.  The  quantity  CT 
should  always  be  the  largest  of  the  three  tensor  contractions 
because  it  is  the  contraction  of  the  full,  nine  component 
gradient  tensor  while  Cl  and  CR  represent  three-component 
subsets  of  the  gradient  tensor.  Differences  in  gradient 
contraction  amplitudes  that  are  measured  by  the  Left  and 
Right  Primary  Guidance  Axes  reflect  the  influence  of  two 
factors  that  affect  the  value  of  the  partial  (sum  of  three 
components)  gradient  contraction  that  is  measured  by  a  given 
sensor  axis.  Namely,  the  magnitude  of  the  partial  gradient 
contraction  C  that  is  measured  by  a  given  axis  (i)  is  a  robust 
function  of  the  axis’s  distance  from  the  target  and  the  length 
of  the  axis’s  projection  in  the  target  direction.  Partial  gradient 
contractions  are  not  full  rotational  invariants;  however,  the 
sensor  geometry  shown  in  Figs.  1  and  2  exploits  the 
symmetry  properties  of  Q  to  develop  robust  magnetic 
anomaly  parameters  that  can  be  used  to  guide  an  UBV  to 
magnetic  targets.  On  the  other  hand,  the  tensor  contraction, 
Ct,  of  the  full,  nine  component  gradient  is  a  true  rotational 
invariant  that  can  provide  an  UBV  with  a  more  robust  and 
longer  range  target-proximity-sensing  capability. 

In  general,  the  scalar  total  field,  BT  (2.2)  is  not  a  robust 
quantity;  however,  as  indicated  by  the  data  of  Figs.  4-6  and  in 
[2],  the  absolute  values  of  the  displacements  of  BTL  and  BTR, 
the  total  field  sums  from,  respectively,  the  Left  and  Right 
Axes  actually  may  be  used  to  provide  a  “pseudo-robust” 
modality  for  magnetic  anomaly  guidance.  The  data  show  that. 


Fig.  3.  Sketch  of  test  range  for  gradiometer  field  tests.  A  permanent  magnet  dipole  target  with  moment  I M  [  =  14.5 
Am2  and  with  declination  of  0°  and  inclination  19°  was  fixed  at  the  center  (0,0)  of  the  17m  diameter  range.  Data 
were  taken  by  moving  the  sensor  along  the  paths  indicated  by  arrows  in  (a)  -(c)-  Part  (a)  represents  “radial” 
approaches  to  target  with  the  sensor  path  passing  through  the  target  position,  while  (b)  shows  “tangential”  approaches 
to  within  1.5  m  from  the  target.  Part  (c)  represents  the  two-dimensional  closure  of  the  tangential  paths  around 
target.  For  the  eight  tangential  paths  of  (b),(c)  the  sensor  system’ s  Left  Guidance  Axis  was  closer  to  target  than  the 
Right  Axis.  Within  the  gradiometer  detection  range,  the  magnitude  of  Left-Axis  gradient  contraction  C  Lwas  larger 
than  that  of  the  Right  Axis’s  Cr  for  all  eight  paths. 
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Fig.  4.  Gradient  contraction,  total  field  and  temperature  data  for  the  sensor  system  on  a  tangential  approach 
to  the  target  with  a  heading  of  135°  (path  6  of  Fig.  3(c)).  For  these  data  curves,  the  left  axis  of  the  sensor 
system  was  closer  to  target  than  the  right  axis. 

The  Total,  Left  and  Right  curves  of  Part  (a)  respectively  refer  to  the  total  gradient  contraction,  CT,  the 
partial  gradient  contraction  for  the  left  axis,  CL,  and  the  partial  gradient  contraction  for  the  right  axis,  CR. 
Similarly,  the  Total,  Left  and  Right  curves  of  Parts  (b)  and  (d)  respectively  refer  to  the  sum  of  total  fields 
measured  by  all  three  TMAs  (BTT  =  BT1  +BT2  +BT3),  the  total  field  sum  from  the  left  axis  TMAs  (BT[  = 
BT1+BT3)  and  the  total  field  sum  from  the  right  axis  TMA  (Bxr  =Bt1  +BT2).  Part  (d)  represents  the  total 
field  data  of  Part(b)  on  a  logarithmic  scale.  For  this  heading  BTL  >  BTR  to  a  distance  of  3.5-4  m. 

Within  the  sensor  detection  range,  the  left  axis  gradient  contraction,  CL,  and  total  field  sums  B-r[  are 
greater  than  the  respective  CR  and  BTR  values  for  the  right  axis.  For  UBV  guidance  to  a  magnetic  target,  the 
condition  CL  >  CR  provides  a  robust  indication  that  the  UBV  should  turn  to  the  left  in  order  to  head  directly 
toward  the  target.  Likewise,  the  condition  BTL  >  BTR  can  provide  a  longer  range  albeit  less  robust  indication 
that  the  target  lies  to  the  left  of  the  UBV  heading. 

Variations  in  gradient  contraction  values  at  distances  greater  than  3  m  are  due  to  TMA  sensor  output 
drift.  The  sensor  drift  appears  to  correlate  with  the  2-3  °C  changes  in  TMA  temperature,  that,  as  shown  in 
Part  (c)  occurred  during  the  data  acquisition  process. 


with  the  Left  Axis  closer  to  target  than  the  Right  Axis, 

Btl  |  >  |  Btr  |  at  target  distances  of  up  to  4  m.  Thus,  for 
some  angles  of  approach  to  target,  pseudo-robust  total  field 
sums  can  provide  longer-range  guidance  and  target  detection 
than  the  robust  gradient  contraction  data. 

As  indicated  in  Fig.  4  and  explicitly  shown  in  Fig.  7, 
sensor  instabilities  and  target-detection-range-reducing  output 
drift  were  found  to  correlate  rather  well  with  variations  in 
TMA  sensor  temperature.  Consequently,  by  solving  the 
temperature-related  sensor  instability  problem,  the 
gradiometer’ s  target  detection  range  and  magnetic  anomaly 
guidance  range  can  be  enhanced.  For  example,  as  shown  in 
Fig.  8,  the  effective  magnetic  target  detection  and  guidance 


ranges  were  effectively  doubled  when  the  sensor  system 
temperatures  were  kept  relatively  constant  (  +  0.5  °C). 

Fig.  9  shows  gradient  contraction  data  that  were  taken 
along  “radial  paths”  such  as  those  indicated  in  Fig.  3  (a). 
These  data  (and  the  data  of  Figs.  4,5,6  &  8)  illustrate  the 
robustness  and  simplicity  of  the  gradient  contraction 
approach  for  magnetic  anomaly  guidance.  For  the  data  sets 
shown  in  Fig  9  (a),  the  sensor  was  pointed  to  the  left  of  the 
target  direction  so  that  the  Right  Axis  closer  to  target  than  the 
Left  Axis  and.  consequently,  CR  >  Cl-  The  condition  CR  > 
Cl  indicates  that  the  sensor  must  turn  to  the  right  in  order  to 
head  toward  the  target.  Fig.  9(b)  presents  data  for  the  same 
radial  paths  as  Fig.  9(a)  but  with  the  sensor 
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Fig.  5.  Gradient  contraction  and  total  field  data  for  a  tangential  approach  to  the  target  on  a  heading  of  180° 
(path  5  of  Fig.  3(c)).  The  Left  Axis  of  the  sensor  system  was  closer  to  target  than  the  Right  Axis  and  the 
magnitudes  of  the  Left  Axis  gradient  contraction  and  total  field  sums  are  greater  than  the  respective  values 
for  the  Right  Axis.  Sensor  gradient  instabilities  were  ss  200-400  nT/m.  Variations  in  TMA  temperatures 
were  very  similar  to  those  of  Fig. 4(c). 


Fig.  6.  Gradient  contraction  and  total  field  data  for  a  tangential  approach  to  the  target  on  a  heading  of  270° 
(path  3  of  Fig.  3(c).  The  Left  axis  of  the  sensor  system  was  closer  to  target  than  the  Right  axis  and  the 
magnitudes  of  the  Left  axis  gradient  contraction  and  total  field  sums  are  greater  than  the  respective  values 
for  the  Right  axis. 


(a)  (b) 

Fig.  7.  Correlation  of  TMA  temperature  and  gradient  contraction  output  drift.  The  gradiometer  sensor  was 
fixed  on  a  heading  of  0°  with  no  target  present  and  TMA-2  of  the  base  axis  was  cooled  for  about  7  minutes. 
Part(a)  shows  the  TMA  temperature  and  Part  (b)  the  corresponding  variation  in  the  Base  axis  gradient 
contraction,  Cb. 
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Fig.  8.  Gradient  contractions  (a)  and  TMA  temperatures  (b)  for  the  gradiometer  on  a  tangential  path  toward 
the  target  located  1  m  down  and  0.28  m  to  the  left  of  the  sensor  path.  The  experimental  conditions  were 
adjusted  to  minimize  temperature  changes.  Consequently,  there  was  relatively  little  drift  in  gradient 
contraction  values  and  the  detection  range  was  enhanced  to  better  than  3  m. 


(a)  (b) 

Fig.  9.  Gradient  contractions  as  the  sensor  system  moved  directly  toward  the  target  (to  a  CPA  of  0.6  m)  on 
radial  paths  of  0°  and  180°  (indicated  by  the  vertical  arrows  in  Fig.  3  (a)). 

In  Part  (a)  the  sensor  was  rotated  45°  to  the  left  of  the  target  bearing  and  so,  at  all  points  along  the  path, 
the  Right  Axis  was  closer  to  target  than  the  Left  Axis.  Consequently,  at  all  points  along  the  radial  paths  of 
Part  (a),  the  gradient  contraction  for  the  Right  Guidance  Axis,  CR,  was  greater  than  the  gradient  contraction 
for  the  Left  guidance  axis.  CL  For  guidance  of  an  UBV  to  a  magnetic  target,  the  condition  CR  >  CL  would 
constitute  an  error  signal  that  indicates  that  the  UBV  must  turn  to  the  right  in  order  to  head  directly  toward 
the  target.  (Conversely,  the  condition  CL  >  CR  would  indicate  that  the  UBV  must  turn  to  the  left  in  order  to 
head  directly  toward  the  target.) 

In  Part  (b)  the  sensor  platform  has  been  rotated  to  point  directly  at  the  target  so  that  the  Left  and  Right 
guidance  axes  are  equidistant  from  the  target.  Consequently  the  gradient  contraction  amplitudes  from  the 
two  axes  are  essentially  equal,  that  is  CR  =  CL,  at  all  points  along  the  paths. 


pointed  directly  at  the  target.  In  Fig.  9  (b)  the  condition  of 
directly  “homing  in  the  target”  is  manifested  by  the  following 
set  of  guidance  parameter  relations:  CR  =  Cl  and  Ct  ,  CR  and 
CL  all  continually  increase  as  the  target  is  approached. 

In  summary,  the  data  of  Figs.  4-9  illustrate  the  following 
characteristics  of  the  sensor  system  of  Fig.  1  and  its 
implementation  of  the  Magnetic  STAR  approach  to  magnetic 
anomaly  guidance  for  UBVs  to  bottom  mines: 


1.  The  gradient  contraction  signal  developed  by  each 
Primary  Guidance  Axis  (i.e.,  CR  and  Cl  )  is  a  robust 
function  of  the  distance  from  each  axis  to  the  target. 
Therefore,  in  order  to  head  directly  toward  the  target,  the 
sensor  platform/UBV  should  turn  in  the  direction  of  the 
axis  with  the  largest  signal  until  the  two  guidance  signal 
amplitudes  are  equal.  As  the  UBV  proceeds  toward  the 
target,  all  C—based  guidance  parameters  should  increase. 
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2.  The  total  gradient  contraction,  Ct,  of  the  full  nine 
component  tensor  will  provide  longer-range  detection 
and  rotationally  invariant  proximity  sensing  for  high 
mobility  sensing  platforms. 

3.  The  sum  of  total  fields  from  each  Primary  Guidance  Axis 
(i.e.,  Btr  and  BTL)  is  a  “pseudo-robust”  function  of  the 
distance  from  each  Axis  to  the  target.  The  Axis  closest  to 
target  usually  will  develop  a  larger  absolute  value  of  total 
field  offset  than  the  furthest  Axis.  This  feature  can  be 
used  to  extend  the  target  detection  and  guidance  range 
beyond  the  gradient  contraction’s  range  [2], 

4.  Sensor  system  instabilities  and  signal  level  drift  appear 
to  strongly  correlate  with  temperature  changes  in  the 
triaxial  magnetometer  field  sensing  elements.  Solution  of 
the  sensor  instability  problem  would  greatly  enhance  the 
effectiveness  of  the  magnetic  anomaly  guidance 
approach. 

The  basic  simplicity  of  the  Magnetic  STAR  approach  to 
magnetic  anomaly  guidance  for  UBVs  is  illustrated  below  in 
Fig.  10. 


Fig.  10.  Simplified  flow  diagram  for  an  algorithm  that  would 
use  signals  from  Primary  Guidance  Axes  and  the  Magnetic 
STAR  method  for  fully  autonomous  guidance  of  UBVs  to 
magnetic  targets. 

V.  CONCLUSION 

The  magnetic  anomaly  guidance  system  for  UBVs  that  is 
described  in  this  paper  represents  an  adaptation  of  the  more 
general,  three-dimensional.  Scalar  Triangulation  and  Ranging 
(STAR)  approach  to  detection,  localization  and  classification 
of  magnetic  targets  such  as  underwater  and  buried  mines. 


Field  test  data  and  a  consideration  of  the  physical  nature  of 
magnetic  anomaly  fields  indicate  that  a  simple,  three-axis 
gradiometer  configuration  will  provide  very  effective,  robust 
and  practical  gradient-contraction-based  magnetic  anomaly 
guidance  for  UBVs  searching  for  mines  on  and  buried  under 
the  sea  bottom. 

Furthermore,  the  power  of  the  STAR  approach  for  DLC  of 
mines  can  be  enhanced  by  increasing  the  number  of  field 
sensing  elements  of  the  array.  For  example,  a  five-primary- 
axis  sensor  system  geometry  can  provide  full,  three- 
dimensional  target  detection,  localization  and  classification 
capabilities.  Such  a  design  could  be  made  sufficiently 
compact  to  provide  a  practical  and  effective  mine  DLC 
system  for  Navy  Divers. 

The  relatively  large  size  of  the  present  sensor  prototype  is 
a  consequence  of  the  instabilities  of  the  “TMA”  field  sensing 
devices.  More  stable  field-sensing  elements,  in  particular 
TMA  devices  that  do  not  exhibit  the  temperature -dependent 
instabilities  that  are  manifest  in  the  field-test  data  would 
allow  development  of  much  more  compact  and  effective 
magnetic  anomaly  guidance  systems  for  high-mobility 
autonomous  sensing  platforms  such  as  UBVs  or  Navy  Divers. 
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